Contemporary Mathematics 



Statistical ensembles and density of states 

Vadim Kostrykin and Robert Schrader 

Dedicated to Yakov Sinai and David Ruelle on the occasion of their 65th birthday 

ABSTRACT. We propose a definition of microcanonical and canonical statistical ensem- 
bles based on the concept of density of states. This definition applies both to the classical 
and the quantum case. For the microcanonical case this allows for a definition of a temper- 
ature and its fluctuation, which might be useful in the theory of mesoscopic systems. In the 
quantum case the concept of density of states applies to one-particle Schrodinger operators, 
in particular to operators with a periodic potential or to random Anderson type models. In 
the case of periodic potentials we show that for the resulting fi-particle system the density 
of states is [(n — i)!2] times differentiable, such that like for classical microcanonical en- 
sembles a (positive) temperature may be defined whenever n > 5. We expect that a similar 
result should also hold for Anderson type models. We also provide the first terms in as- 
ymptotic expansions of thermodynamic quantities at large energies for the microcanonical 
ensemble and at large temperatures for the canonical ensemble. A comparison shows that 
then both formulations asymptotically give the same results. 



1. Introduction 

The classical concepts of canonical and grand canonical statistical ensembles are well 
known to have important quantum analogs (see, e.g., standard textbooks like [^p[|35], |64|]). 
To the best of our knowledge the notion of a microcanonical ensemble in quantum theory 
has not received the attention we think is deserves. This is unfortunate in view of recent 
developments in the classical and quantum theory of small ensembles. So far the main mo- 
tivation for a microcanonical analysis of small, classical systems came from gravitational 
physics with its long range forces (see, e.g., ^, Bl, 75, 3^, |52|, ^] and the references 



quoted there) and from nuclear physics with the associated fragmentation processes (see, 
e.g., [ |28| , [2^] and the references therein). The main aim was to ask for the most probable 
distribution and to view this as an equilibrium phenomenon. By energy conservation for 
closed systems it is natural to consider equilibria in the microcanonical description, where 
the temperature T and its fluctuation are derived quantities. For the operational definition 
of a temperature for small systems one cannot use a heat bath but rather a small thermome- 
ter, see [7|, p.351]. 
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Usually one introduces the notion of statistical ensembles to describe the systems with 
large number of particles. Our approach will in some sense go in a different direction. We 
deal with the question how small the system can be to allow a thermodynamic description. 

We recall that in the familiar classical microcanonical description IIT may be written 
as a mean over the energy shell of a certain function in the (finite dimensional) phase space 
(see relation (2.5) below). Considering the square of this function one may in addition 
consider fluctuations of the inverse temperature. In case the ergodic hypothesis holds such 
averages may then be written as time averages using the dynamics of the classical Hamilton 
function (see, e.g., [ |67| , |68[ |L ^|). An important feature of a microcanonical ensemble is 
that the temperature may decrease with energy giving rise to a negative heat capacity. The 
first references which envisage such a situation seem to be [^0| and pfl|]. The example 
considered in [ SO | is that of a star, which has used up its nuclear fuel and then heats up 
under gravitational contraction. A first model discussion for a supernovae providing a 
negative specific heat was given in [^5[. For recent experiments on small systems, see, 
e.g., [|28[ p9|, [7(^ and references given there. For more recent theoretical discussions of 
small systems, see, e.g., $i% ^4| . We mention also the recent work [^4[ where a system 
of classical particles with power law potentials was considered. The results of this work 
indicate that already for particle numbers larger than 10 there might be a relaxation to the 
microcanonical equilibrium starting from an arbitrary initial state. 

It is the purpose of this article to provide a useful quantum theoretical approach to the 
theory of canonical and microcanonical ensembles which may be applied to small (i.e., 
mesoscopic) systems. This will include a discussion of fluctuations, the heat capacity 
and their interrelation. We also compare the microcanonical and the canonical ensemble. 
Our approach will be based on the theory of (integrated) density of states, which in the 
quantum case is a well known notion in solid state physics (see, e.g., [Q, |55|]). This theory 
has been investigated intensively in the last years by mathematical physicists in the context 
of one-particle Schrodinger operators, in particular, operators with a periodic potential and 
random Schrodinger operators like Anderson type models. For recent references consult, 
e.g., [gggg] and for references before 1992 [g^^^]. 

In our approach the microcanonical description is given by defining the entropy to be 
k (the Boltzmann constant) times the logarithm of the density of states, so that the entropy 
comes to be a function of the energy. As usual the inverse temperature IIT is then the 
derivative of the entropy, i.e., k times the logarithmic derivative of the density of states. 
The canonical description on the other hand is by definition given in terms of a partition 
function, now defined to be the Stieltjes-Laplace transform of the integrated density of 
states w.r.t. the variable p = llkT canonically conjugate to the energy. 

In particular we will address the familiar question under which conditions these two 
descriptions give approximately the same answer. In fact, we will see that the well known 
methods of comparison in the standard formulation easily carry over to this new formu- 
lation. In addition we will show, both in the classical and quantum cases, that for high 
energies or correspondingly high temperatures the microcanonical and the canonical de- 
scription give the same results. 

So for our approach to work in the microcanonical case, differentiability of the inte- 
grated density of states up to third order is necessary. Now our main observation is that 
smoothness of the integrated density of states for the resulting n-particle theory increases 
with n. This will in particular allow us to consider microcanonical ensembles for particles 
moving in a periodic external potential provided n>5. The reason is that for one-particle 
Schrodinger operators with a periodic potential the density of states (which is the derivative 
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of the integrated density of states) is well known to exhibit so called van Hove singulari- 
ties in space dimensions greater or equal to two (see, e.g., [|[ |55|]). In one dimension the 
integrated density of states itself has square root singularities. These are smoothed out 
when going to higher particle number. We expect a similar property to hold for particles 
moving in an external random potential. We will argue why increased smoothness with 
n is related to the well known increasing smoothness of the integrated density of states 
with the dimension of the space. For Anderson-type random Schrodinger operators the 
existence of the density of states (i.e., the absolute continuity of the integrated density of 
states) is a complicated, still not completely solved problem. Very little known is known 
about regularity properties of the density of states [ IjjL p3| , BgLBjh. 

In addition we will invoke notions from scattering theory and, in particular, the scat- 
tering phase shift (or total scattering phase) at a fixed energy E. It is an old observation of 
Beth and Uhlenbeck [|7| that the second virial coefficient i n stat istical mechanics is related 
to the phase shift (see also [[35]] and relations ( |3lj|) and ( J3.24 ) below). We will use the 
concept of the phase shift density introduced in [|45|, [4(j | to give the asymptotic behavior 
at large energies (or temperatures) of the main thermodynamic quantities. In particular, in 
[ fis] , [46[ | we established that up to a factor of 7t this phase shift density is the difference of 
the integrated density of states with potential and the free theory. This relation is analo- 
gous to the change of the number of particle states found by Friedel [ |20[ ^T] l in the case of 
a single impurity. The phase shift density will be used to establish relations on the shift of 
the temperature induced by the potential in the microcanonical context and on the shift of 
the mean energy in the canonical context. Furthermore, we establish a relation on the shift 
of the mean energy density for a system of noninteracting electrons moving in a periodic 
or a random potential. This relation is analogous to a theorem of Fumi [ |23| ] which relates 
the shift of the ground state energy due to a single impurity to the scattering phase. 

We will not discuss situations where in addition to energy other quantities, like, e.g., 
angular momentum, are conserved. Also we will not cover the situation where in addition 
to an external potential there is interaction between the particles. This is of course an 
important issue worth pursuing. We note that non-relativistic n-particle scattering theory 
with all possible fragmentation and bound state channels has been extensively analyzed 
(see, e.g., [ [ffi p^ , |36| ] and the references quoted there). We expect that some of our results 
may be extended to this situation, although at present it is unclear to us, which role such 
multichannel spectral properties play in the theory of integrated density of states. 

In addition, it would be interesting to see whether the present approach to the theory 
of quantum mechanical, microcanonical ensembles could serve as a laboratory for a fresh 
look at ergodic theory in quantum mechanics (see, e.g., [^8, 62, 19|). 

The article is organized as follows. In Section g we will review the classical theory 
of microcanonical and canonical ensembles in a form which serves as a motivation and 
comparison for our quantum mechanical approach to be given in Section |[ For the clas- 
sical models we consider, we will introduce a classical notion of (integrated) density of 
states and formulate the resulting theory of statistical ensembles. We will see that in this 
classical setup differentiability indeed increases with the particle number. In addition, we 
will provide a new example for a system with negative heat capacity. Also we will provide 
an example where the microcanonical and the canonical descriptions asymptotically for 
large E and T give the same (mean) energy - temperature relation. Recall that in the usual 
canonical description the heat capacity is always positive. 

So the discussion in Section |2] will be rather extensive, the reason being that the quan- 
tum theory is then simply obtained by replacing the classical (integrated) density of states 
by the corresponding quantum version. 
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Section [3| will contain the corresponding quantum mechanical formulation of both the 
microcanonical and the canonical ensemble using the quantum notion of the integrated 
density of states. For large energies we again will argue on the basis of a yet unproven 
conjecture on the high energy behavior of the phase shift [ f7l| ] that the (mean) energy - 
temperature relations are asymptotically equal for both ensembles. As already mentioned 
we view this result as an indication of the reasonableness of our thermodynamic approach 
to small systems. Section || will also include a brief comparison between the classical and 
the quantum theory for small h. Finally in Section |3| we also briefly discuss the grand 
canonical ensemble in terms of the integrated density of states. 

In the appendix we will show in the quantum context for periodic potentials that sin- 
gularities of the n-particle density of states are smoothed out with increasing number of 
particles n. In addition we will show in the classical case how the randomness in stochastic 
potentials smoothes out the integrated density of states. 

We avoid to formulate our results as theorems and propositions. However making any 
statement we provide a discussion of what is rigorously proven and what is only conjec- 
tured motivated by the physical intuition. 

Acknowledgments. We are indebted to R.M. Dudley, P.D. Hislop, M. Karowski, 
H. Narnhofer, M. Schmidt, K.-D. Schotte, H. Spohn, and W. Thirring for valuable com- 
ments. We thank D.H.E. Gross and O. Fliegans for providing references. 

2. The classical theory 

In this section and for the purpose of comparison and as part of our motivation we 
will briefly review the well known concepts of microcanonical and canonical ensembles in 
the classical case. Some of the material, however, seems to be new. In particular we will 
provide new examples with negative heat capacity. Previous examples in the context of 
gravitation theory are given, e.g., []52|,|75[] and the references quoted there. 

In classical theory the starting point is a finite dimensional phase space 7 which is 
a symplectic manifold. We will assume, in addition, that 7 also is a Kahler manifold, 
i.e., it is also a Riemannian manifold and the symplectic and Riemanian structures are 
compatible. This implies, in particular, that the Liouville volume form and Riemannian 
volume form agree. We denote this volume form by livoly. Also a Hamilton function H c \ 
on 7 is supposed to be smooth and bounded below in case 7 is not compact. 

First we assume f<E = {/)£ CP | H c \{p) < E} is compact for all E. Note that Hamilton 
functions of the form// c i(p, x)=p 2 /2m+V(x) with bounded V do not satisfy this assumption. 
We will turn to this later. 

The function 

(2.1) N(E;H c i)= [ dvoly = [ @(E -H^dvoly, 

Jy< E J? 

where © is the Heaviside step function, is increasing w.r.t. the energy E. Its derivative is 
given as 

< W(E;H cl ) = —N(E; H cl ) = [ 8(E - H cl )dvol? 
dE Jy 

with 8 being the Dirac 8-function. The associated microcanonical ensemble is then given 
by defining the entropy as a function of the energy E as S(E,H C \) = kin W(E;H C {) and 
the temperature T(E) = T(E;H C {), again considered as a function of the energy E, is then 
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defined by 

1 dS(E;H cl ) d dN(E;H cl ) 1 dW(E;H cl ) 

(2.2) = =k — In = k 

T(E) dE dE dE W(E;H c i) dE 



with k being the Boltzmann constant. So provided the r.h.s. of ( |2.2| ) is meaningful the tem- 
perature is defined but possibly negative. In standard situations W(E; H c \) increases with 
E such that T is nonnegative. However, there are also situations, where the temperature T 
may become negative (see, e.g., [pfl]). 



The function W(E; H c {) and its two first derivatives may be written in the form 

W(E;H cl ) = f dfi E 

d f V# c i 

(2.3) ^W(E;H cl ) = / V • — -^dp E 

dE Jy E IVHcil 7 

with the following notation. The manifold CPg, the boundary of 7<e, is the energy surface 
for the energy E, i.e., the set of points p e V for which H c \(p) = E. The measure dps on 
Ve is given as dps = \VH c \\~ l dvoly E , where d\oly E is the canonical volume form on Te 
given by the Riemannian metric. V is the covariant gradient and |Vi/ c i|(/?) the length of 
the vector VH c \(p), which is normal to the surface Cp£ at p e y E . So provided VH c i(p) # 0, 
the vector l/\VH c i(p)\ VH c \(p) is the outward unit normal vector to 7 E at p e 3>e, i.e., 
this vector points into the complement 7>e of IP<£. Near Te one has the familiar relation 



afvoly = dp.EdE, whence the first relation in (2.3). Finally V- is the covariant divergence. 



The second and third relation in (2.3) are obtained from the first by (repeated) use of the 
theorem of Gauss (see (2.11) below). Higher derivatives may be calculated similarly. In 
standard situations is a compact set, generically of codimension 1. For W c i(E; H c {) and 
its derivatives to be well defined one has to assume that VH c \(p) £ for p e Ve (or at least 
that suitable inverse powers of \VH c \\(p) are integrable over Te)- Then in particular is 
a smooth submanifold of V of codimension 1 and W(E; H c i) ? provided the energy E is 
such that y E * 0. The first relation in (O) dates back to Khinchin [B2|]. 



The temperature T(E\H C \) defined by (2.2) is an increasing function at E if 



(2.4) W(E;H cl )-^W(E;H cl )<(-^W(E;H cl ) 

We may rephrase this as follows (see also |]67|, |6^]). Let 

, , _ U E fdVE _ J 5(£ - H ci )fdvol v 
{f)E Jj, E dp E J8(E-H d )dvoh 

denote the average over the energy shell y E of a real valued function / defined on y E . For 
given H c \ consider in particular the function T on 3 of the dimension of energy given as 

r' = y. VHcl 
|VHaP 

such that 

(2.5) m =k{7 ~ l)E - 

In general kT(E) ^ (T)e- Consider for example the case where < T< °o on y E , Then 

\ = {\) E = {7 m T m )E<{T ) f(T l )f 
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by the Schwarz inequality giving (7)e > kT(E) with equality if and only if T is constant 
onT £ . 



We write 



(2.6) 



V V 



V// C l 



|VH.j| : 



1 

T 2 



iVfldl 2 



+ s. 



: (VT 1 ) 



Here • denotes the scalar product of two vectors (or rather vector fields) given by the 
Riemannian metric, such that in particular |Vi/ c i| 2 = VH C \ ■ VH c i. 
Consider the quantity A(r _1 )(£) = A(T~ l )(E;H c \) given by 



(2.7) 



l )(E) 2 -- 



:k 2 ((T 2 ) E -(T l ) 2 E )>0 



may be viewed as the fluctuation of the inverse temperature at energy E. Note that the in- 
equality is a consequence of the Schwarz inequality or equivalently of the familiar relation 
< ((/ - (f)E) 2 )E- In particular A(T~ l )(E) vanishes if and only if T is constant on CPg. 
The fluctuation and the heat capacity c v (E) = c v (E; H c \) defined by 



(2 - 8) - tW 
are related by 

(2.9) 



dE 



r, d 1 r, d d 

T{E) = -T(E) 2 —— = -kT(E) — In — W(E; H cl ) 
dE T(E) dE dE 



A(T- 1 )(E) 2 + 



T(E) 2 c v (E) 



-k 2 {5)t 



Therefore condition (2.4), which guarantees that c v (E) > 0, may be recast into the equiva- 
lent form 



(2.10) 



- l )(Ef < 



■<S)i 



By we have \VH ci \5 = V(l/T) • |Vi/ c i| _I VH ch This has the following geometric inter- 
pretation. | V// c i|S is the component of the gradient of 1/T in the normal direction. So if 
this gradient always points into !P<£ (i.e., V(l/T) • VH c \(s) < for all s e Te) thus making 
9 negative (or zero) there, then the inequality (2.1C) is satisfied for the energy E. Since 
V(l/T) = -(1/T) 2 V T this is equivalent to the conditio n that the gradient of T always points 
into 7 >E = y\ y<E- Conversely, if the inequality in ( 2.1C ) is reversed and in particular if 
(5)e > (meaning that in the mean the gradient of 1/T points into O 3 ^, then T(E, H c \) is a 
monotone decreasing function in E, i.e., the heat c apac ity is negative. 

By means of the Gauss theorem and using (2^) the relations (23) can be written 
equivalently as follows 



(2.11) 



W(E;H d ) 



—W(E;H ci ) 
dE 



T l d\oh, 



(T l +7Q) T l dvoh, 



dE 2 



W(E;H d ) 



T" 2 + 2g + Tvg • 



T-'Jvoly. 



Moreover, for any differentiable function / on T<£ the average over the energy shell can 
be written as 



(2.12) 



;/)/ 



+fr 1 ) dvoh 



J 9 J- l d\olj> 
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Assume now that 7 > on 7<e- Then we may introduce the probability measure on 

3><e by 



(2.13) 



1 



W(E;H d ) 



©(E-H^T^voly. 



In analogy to {-)e let {-)' E denote the resulting mean. Then relations (2^) can be rewritten 

as 



1 



■■k{T l +7S)' E . 



T(E; H c \) 

However, the resulting inverse temperature fluctuation given by 

-h/i7i _ /;,2//<r-l , <tpi2\/ tv j7i-2\ 1/2 



(2.14) 



A'(r-')(£) = (k 2 ((7~ l + 75) 2 )' E - T(E)~ 



is in general different from A(T l )(E). In fact, by (2.12) and the Gauss theorem 

< (T 2 } E = ((7~ 2 + 23)}^ = ((7~ l + 7Q) 2 )' E - ((75) 2 )' E . 

This implies A(T )(E) < A'(T~ 1 )(E) which is an equality if and only if TS vanishes 
identically on CP<£. Relation (2.9) is now replaced by 



A'(T- l )(E) 2 + 



T(E) 2 c v (E) 



■ (:rvs -&' 



so alternatively to ( J2. 1 0| ) the condition ( |2.4| ) for the heat capacity to be positive can now be 
written as 

I V/-/ 



The alternative definition (2.14) of the inverse temperature fluctuation has the disad- 
vantage that besides the assumed positivity of T it involves the values of H c \ on all of CP<£. 
In contrast and due to the appearance of derivatives of H c \, the definition (2.7) only involves 
the values of H c \ near Cp£. 

The Boltzmann-Gibbs partition function in the associated canonical ensemble may 
also be given in terms of the functions N(E; H c i) and W(E; H c i) as 

(2.15) 

Z(p;// Cl ) = J e~ fiHd dvoh = j e-V E dN(E;H cl ) = j e^ E W(E; H cX )dE, $=l/kT. 
The (mean) energy-temperature relation now takes the form 



£(p; ff cl ) = -^Z(P; H c i) = (£}p with 



Jg(E)e-V E W(E;Hd)dE 
Z(p;// Cl ) 



being the mean of any function g of the energy for given p. The mean energy now is a 
monotone increasing function of temperature. Since (E 2 )$ > {E)l by the Schwarz inequal- 
ity, A£(p; H c \) = ({E 2 )p - (E)p) 112 gives the energy fluctuation. Stated differently, the heat 
capacity defined in this canonical context as 



(2.16) 

is always nonnegative. 



cv(p; ffd) = ^£(P; tfd) = ~ ^£(P; H cl ) 
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The following discussion (see also [68|) relates the microcanonical and canonical de- 
scriptions. It is an adaptation of well known arguments employed in this context (see, e.g., 
[Ba|) and is recalled for later purpose. Write 



Z(P;ffd)= ( e-^ w(E ^dE. 



(2.17) 



The integrand takes an extremal value at E m . dx = E max (T) given implicitly by 



(2.18) 



1 

T 



k— \nW(E;H cl ) 
dE 



E—E m 



which agre es w ith ( |2.2| ). If a solution £ mal to the relation ( |2.18| ) exists it is unique. If the 



inequality ( |2.4[ ) holds for all E and is a strict inequality at E = £ mal then the extremal value 



of the integrand is a maximal one. Indeed, for E near E max we obtain 



1 



-p£ + ln ( ^W(E:H, } ) 



1 



-ME - E max f + 0((E - £ max ) 3 ) 



with a = a(p; //<"') given by a = -d 2 In W(E = £ max ; H d )/dE 2 , i.e., 



(2.19) 



dE 



■W(E;H cl )-W(E;Hci) 



W(E;Hc\) 2 



which is nonnegative if (2.4) holds and is strictly positive if (2.4) is a strict inequality 
at E = Emax. In fact, in the m icroca nonical notation a is related to the heat capacity via 
a= ll(kc v (E m . dx )T (E mlix ) 2 ) by (2.18). So when a > is large, to a good approximation the 



integrand in ( 2.17 ) is given by a Gaussian distribution with variance a . The energy flue 



tuation is then A£(p;H c i) ~ a~ I/2 . Finally if in the canonical ensemble the heat capacity 
is defined as c v (T) = dE($; H c \)ldT this quantity then agrees with the microcanonical heat 
capacity at the energy E mdx , i.e., kT 2 c^(T) = A£(p;// C i) 2 ~ kT 2 c v (E mdx ). 

^From now on we will consider classical Hamilton functions with non-compact T<e 
such that the above discussion does nor apply. The following discussion serves as a prepa- 
ration to what we will do in the quantum case in the next section. More concretely, we 
consider a particle of mass m > moving in R. d under the influence of a periodic bounded 
potential V(x), i.e., V(x) = V(x+ j) for allxe W l and all je Z d . Thus the classical Hamilton 
function is given as H c \(x, p) = p 2 l2m + V (x) and the phase space 7 is M. 2d with the canoni- 
cal symplectic structure. We will consider the resulting n-particle theory with phase space 
M. 2nd and classical Hamilton function H$\x, p) = p 2 /2m + V {n \x) = H^(p) + V {n \x). Here 
p = (p\,p2, ■ ■ ■ ,p»),x= (x\,X2, ■ ■ ■ ,x„), pi.xiG R d and V { "\x) = V (x,). For the purpose 
of this article it will suffice to assume the potential to be bounded. 

We will also consider random, classical Hamiltonians of the so called Anderson type. 
More precisely, we assume the one particle random potential V = to be of the form 



(2.20) 



V<aW= £ qj((0)Vo(x-j)- 



Here C0={<D j , } j(e z<' is an element of the probability space £2=]R Z with a probability measure 
Prob and qj((0) = q((Oj) for some real valued bounded measurable function q. Therefore 
CO/ are independent random variables with the same distribution. In other words there is a 
probability measure d/j on M such that Prob(^j £/) = /, dp(q) for any interval / cz R and 
any je 1/ . Also we assume that Vb is a bounded function with support in the unit cube in 
M. d , such that Vq(x - j) and Vo(x - k) have non-overlapping support when j + k. This results 
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in a Hamiltonian H c \(p,x) = H c i a (p,x) = p 2 l2m + V®\x), again with the obvious notation 

v^ ) M = L-Va,(* ; ). 



In the considered cases the function N(E;H^) defined by (2.1) is infinite but the 
integrated density of states, given as 

n(£; H ( " ') = lim — - / ®(E - H%\p, x)) d nd x d" d p 
a^» |A|" jA»m m > 

exists where A c R d and A" = A x . . . x A c W' d . If V is periodic we obviously have 
n(£; //<;")= / ®(E-H^(p,x))d nd xd" d p 



(2.21) " 

= —Vol(S" d - [ )(2m)" dl2 / ®(E - V (n \x)){E - V (n \x)) ndl2 d nd x. 
nd JA'i 

where Ao is the unit cube in R d , §" rf_1 is the unit sphere in W d and Vol(S" fl,_1 ) its area. In 
the case of free motion, i.e., V = 0, and using the fact Vo^S" 11 " 1 ) = 2% ndl2 IY(ndl2) we obtain 

(2.22) n(£; ff£>) = — Vol(S" rf - 1 )(2m£)'" ,/2 = { - 2%mE ^ m E > . 
J V ()ci; nd T{ndl2+\) 



Note that from (2.21) it follows that the order of differentiability of n(E;H^) w.r.t. E 
increases with nd. 



In the random case V = V m we obtain in place of (2.21) 
n(E;H l f) = -i ; Vol(§' , ''- I )(2m)" "" 



(2.23) Ud 



[ dp(q) [ ®(E-qT V (Xi))(E -qY V ( Xi )y d > 2 d" d x. 
JR JA"„ i i 



For fixed number of particles n the smoothness of n(E;H^') depends on the probability 
distribution dp(q) of random variables qj(oS). Assume that the measure dp(q) is absolutely 
continuous with smooth density supported in an interval [e, £ _1 ] for some < £ < 1. Then 
we may rewrite (2.23) as 



n(E;H { " ] ) = — Vol(S' tt/ - I )(2m) m ' /2 
nd 



(2.24) 



j q" dl2 f(q)dq J © (j- - V (x)\ V (x)\ d 



nd/2 

■nd t 



j„d/2+i J_Vol(S" c '- 1 )(2m)'" </2 
nd 



( q" d/2 f(Eq)dq ( ( - 
J JA a \q 



I \ nd/2 



q V °W)[-- Vq W) d " dx - 

In particular, if / e C°° then the degree of differentiability of n(E;H^) is completely 
determined by the factor E" d/2+l . 

The main idea of the present work is to replace N(E;H^) in the relations (2.1) - 
(2.2) by the integrated density of states n(E;H^'). For nd > 3 this gives the temperature 
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T(E) = T(E;H% ) ) as 



1 d d i„\ 

— n — n(E;H%') 



T(E) dE dE 

(nd \ $ K ®(E-V {n Xx))(E-V {n) (x)) ndll - 2 d nd x 
(225) = \J~ 1 ) k J^ ®(E - VW(jt))(£ - y(n\ x ))«dii-\ d nd x ' nd ~ 3 - 

which is obviously nonnegative. 

The motivation is as follows. We consider the periodic case only, the random case can 
be treated simi larly . For a given finite volume A the temperature can be defined by means 
of the relation ( |2.2| ), i.e., 

1 k ± ln dN ^ H ^ 



T(E; H^ A ) dE dE 
with 

N(E;H^ A )= [ ®(E-V M (x))(E-V M (x)) ndl2 - 1 d nd x. 
Ja" 

Considering the limit A — » °° in the Fisher sense due to the relations 

4-jn(E;H^)= lim f t N(E;H%\), jfc=l,2, 
dE k cl a^oc |A|» dE k dA 

we obtain lim T(E; H%\) = T(E; H§\ 
The condition 

d ,„ 'T<-"h—n(E-H (n) ) < (— -'^ 
dE— 1 * >dEi Ht ' Hcl) - [d0 



(2.26) — n(£;^) —n(E;H%>)< I n(E;H%')) , nd>5 



guarantees that the temperature T(E; H^') is a non-decreasing function of the energy, i.e., 
it naturally replaces condition (2.4). It may be rewritten as 

(2.27) (^-2)(T £ -) £ ,(^-l)(T £ -')| 

Here T £ (x) _1 = (nd/2- l)(E -V (n) (x)) and ( ) E is the mean given by the probability 
measure 

&(E - V {n) (x))(E - V in) (x)) nd,2 - 1 d nd x 



dfj E (x) ■■ 



J A „ ®(E - V {n \x))(E - VM(x)) ndl2 - 1 d nd x 



on Aq provided E is such th at th e set {jt|V' n) (jic) < E} has non-zero measure. Thus, we 
obtain the relation similar to (2.5), 

(2.28) _1_ = A:(T £ - I } £ . 



Note that the heat capacity will be negative if the inequality opposite to ( |2.27| ) holds. 
By the Schwarz inequality we have 

(2.29) (7 E - 2 ) E >{7 E - 1 ) 2 E . 

On the other hand (nd/2 - 2) < (nd/2 - 1). Therefore for the heat capacity to be negative the 
inequality (2.29) has to be sufficiently strict. 

As an example where the temperature decreases with the energy consider the following 
choice for the periodic potential. Let XcAo have measure < a < 1 . Define V (x) = Eq > 
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for xe X and zero otherwise when x e Aq \ X. Finally extend V periodically to all of WL d . 
In what follows Eq will be fixed. Then for any E > uEq and k = 1,2,3 



0< / &(E-V {n) (x))(E-V {n) (x)) ndl2 - k d nd x 



(2.30) 



E 



a J (l-a)"- J (E-jE } 



nd/2-k 



We now fix a to be given as a = 1 - ((E - nE^)IE~) (nda ~ T,ln . Then the leading contribution 
to (2.30) when E - hEq is small is given as 



(E 

const • (E - 



nEoY 
nE ) 



ndl2-2 



for 
for 



-3, 

: 1,2. 



When inserted into (2.27) the l.h.s. behaves to leading order as const • (E - nEo) ndl2 ~ 5 , 
whereas the r.h.s. behaves like const • (E - nEo)" dl2 ~ 4 with positive constants depending on 
E, n and d. Therefore when E > nEo is sufficiently close to nE the inequality in (2.27) 
is indeed reversed. Having thus fixed a by continuity the heat capacity therefore also 
becomes negative for all energies sufficiently close to E. Starting from the potential just 
constructed by a small change we may achieve that V is smooth. By the same arguments 
we may replace A () by any A and find a potential V supported in A such that the heat 
capacity defined by N(E; ffi A ) is negative at least in some energy interval. 

By the representation ( 2.28 ) we may also define the inverse temperature fluctuation 
A(T- 1 )(E) = A(T- 1 )(E;H^ ) )by 

(2.31) A(T- 1 )(E) 2 = k 2 ((7 E - 2 ) E -(7 E 



-1\2* 



> 0. 



Recalling the definition of the heat capacity as llc v (E) = dT(E)ldE we obtain in the case 
nd>5 



1 



cAE) 



-kT(E) 



( nil 
V l 



i nd 
v 2 



2) (Te- 2 )e- 



1) 



and therefore 



(2.32) 



A(T- 1 )(E) 2 + 



nd 
~1 



i Cr E - 2 ) E >o. 



T(E) 2 c v (E) 

In the absence of external fields, i.e., for V = this gives A(llT)(E) = 0. Since this is an 
undesirable feature, we will consider the following alternative. Let E^- n = inf teR ,„i V (n \x) > 
—oo. Motivated by the construction (2.13) consider the probability measure on the interval 



(«) 



(2.33) 



cN E (E') ■■ 



1 



U(E;Hf)dE 



jn(E';H™)dE' 



dE 



such that all derivatives (provided they exist) of n(E';H^ ) with respect to E' > E„' n are 
positive. Denote the resulting mean by ( )' E . Observe that 

dv E (E') 



dE' 



E'=E-0 



kT(E) 



and so again we have 



T(E) 



-k(T l y E 
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with 



, -i d d 

7(E') 1 = In 

V ; dE' dE> 



-n{E'-Hf), 



In analogy to (2.14) this results in an alternative definition of the fluctuation of the inverse 
temperature of the form 

A'(T- 1 )(E)=(k 2 (7- 2 Y E -T(E)- 2 ) m 



and which in general differs from A(T l )(E) given by (2.31). For comparison we first 
observe that (nd > 5) 



(^-1) d* 



ME\H^) I —n(£; 
(f-2)dE i cl I dE cl 



On the other hand 



(2.34) (T 2 )' E 



^(E'-H^) 



dE 



Tl n(E';H^))dE' I ^-n{E;H^). 



(«K 



dE 



dE 



For nd > 7 we perform a partial integration in (2.34) and obtain by comparison 

(— - 2) 

(235) (T-y F = ( -±>(7 E - 2 ) E+ (S) E 



with 



Cry = ul 

\ IE { nd_ V) 



~S(E')- 



i 



dE' 7(E') 

The condition S(E') > 0, which means that 7(E') increases at E', is equivalent to 



n(E';HM)^n(E';H™) < [ ^n(E';H%>) 



dE' 



dE 



which compares with (2.26) and hence may be discussed in a similar way. By (2.32) and 
( 2.35 ) the heat capacity can be expressed as 

1 



er 2 y E -{9y E -(7-') 



E 



This alternative definition of the fluctuation and the heat capacity are therefore related by 

k 

-' T ' rW "/£ 



(2.36) 



A'(r _1 )(£) 2 + . 



^r(£) 2 (S)i 



T(E) 2 c v (E) 

As a test in the special case V = 0, i.e., in the absence of external fields, and for nd>6 
this new fluctuation is calculated explicitly to be 



(2.37) 



A'iT-'XE)-- 



nd 



-1/2 



nd 



-1/2 



1 



T(E) 



If we define the temperature fluctuation as A'(T)(E) = A'(T~ l )(E) T(E) 2 then for large nd 
we obtain A'(T)(E) ks (2lnd)T(E), now as it should be. 

The following asymptotic expansion holds for large E and nd>2 



(2.38) 



r(-2" + 1) d („) nd 

(2TZm)^ l2 dE n{E ' McX )= ~2 
nd 



®(E - V (n \x))(E - V m {x)T alA - l d na x 



x, 



W I „\Y'd/2-l jnd 



_ E „d/2-l j-j _ aE -l + fe£ -2 + (9(^-3)] 
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where 



(y (n) W) 2 c/" d x 



nV 2 + n(n - 1)V 2 



Here we use the following notation. In the periodic case and for any natural number k 



(2.39) 

In the stochastic case 
(2.40) 



yk, 



V(x) k d d x= / V {xfd d x. 

A J h 



V k = q k / V (xfd d X with q k - 



q dfj(q). 



Note that V~ < V 2 by the Schwarz inequality with equality if and only if V is constant in 
x (and in co in the random case). This gives the following asymptotic energy-temperature 
relation 
(2.41) 



1 



T(E;H™) 



' nd 

-M--UE- 



l+E~ l nV +E~ 2 n 



+ 0(E~ 3 ) 



Therefore for all sufficiently large E and for nd> 2 the temperature increases with E. Also 
when nd> 2 at fixed large energy E to leading order in E~ l the temperature T(E) decreases 
(as a function of the potential) when a potential with V > is switched on. Otherwise it 
decreases. Similarly if V = such that V is not a constant, then T(E) decreases for all large 
nd when such a potential is switched on. 

For the heat capacity we obtain the asymptotic expansion for large energies 



(2.42) 



c{E;H^)-- 



nd 
~2 



1+E~ 



nd 
~2 



v 2 -T 



+ 0(£" J ) 



Thus for large E and nd>4 the heat capacity increases when a potential is switched on. 

A similar discussion for a canonical ensemble is possible by defining the partition 
function as the Laplace-Stieltjes transform of the integrated density of states 



e-V E dn(E;HW) = $ / e^ E n(E;H^)dE 



such that z(p; Hq'^) = (2%mkT)" dl2 . This limit has not to be confused with the notion of 
"thermodynamic limit" used when discussing extensive systems in statistical mechanics. 
There one divides the logarithm of the partition function in a finite volume by the volume 
of the configuration space and then takes the limit (see, e.g., [^|). So in order to distin- 
guish these two limits we use the notion "per unit interaction volume" (see, however, the 
discussion of the quantum mechanical grand canonical ensemble in Section ||). The reason 
for this difference is as follows. The thermodynamic limit has to be taken when the number 
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of particles increases with A such that the density remains constant. Here the system is not 
extensive since the number of particles stays fixed. 

The resulting entropy, mean energy and the Helmholtz free energy are then given as 

-kf-^dz (ftfl* )), 



In z(p; //<;>), 



/(Prf) 



1 



"P 



lnz(P;// c ( ;°) 



giving the familiar relation between these three quantities 



(2.43) 



e -(P;^)=^(P;< J )+/(P;^" J )- 



In analogy to ( 2.16| ) the heat capacity in this canonical setup is given as c v (T;H^ c f) 
de($;H ( f)/dT. 



Also in analogy to the discussion following (2.17) we may relate the microcanonical 
and the canonical description. We will assume nd > 5 such that n(E,H d ) is three times 
differentiable w.r.t. E. Write 

^^-n(E;H^)dE: 



dE 



(2.44) z(P;#™) 
Then the integrand takes its extremal value at E, 
(2.45) 



P 



-P£+In(i 



7 m\Hf)) dE _ 



' ^max 

(T) given implicitly by 



T dE dE cl ' 



\E=E m 



Relation (2.45) now compares with (2.25). In other words the temperature defined in the 
microcanonical ensemble at the energy E=E max agrees with the chosen one for the canoni- 
cal ensemble. If a solution E mdx to this relation exists it is unique and a maximum provided 
the inequality (2.27) holds for all E (for all large E this is true by (2.42)) and is a strict 
inequality at E = £ max . Indeed, for E near £ max we obtain 



' E mdx ) 0({E E max ) ) 



with a = 



: a(p; ffgi ) now given by 



a = 



on 



E-E m 



which is nonnegative if ( 2.27 ) hol ds an d is positive if ( |2.27| ) is a strict inequality at£'= £' max . 
Again as in the discussion after (2.17) and in the microcanonical notation, by (2.45) a is 
related to the heat capacity via a = l/(Ar v (£ max )r(£ max ) 2 ). 

Hence to a good approximation and for given temperature the canonical ensemble 
is described in terms of a Gaussian distribution in energy with variance a~ 1/2 in case a 
is positive and large. In particular in this approximation we obtain e(P;// c ( "') w E m - dx 
and Ae($;H^) w a~ m . Similarly the canonical heat capacity and the microcanonical 
heat capacity at energy £ mal agree in this approximation (see the discussion after (2.19)). 
Observe that for large T the implicit equa tion (J2.45 ) has always at least one solution E max w 
(nd/2 - \)kT due to (2.38), which by (2.42) is at least a local maximum. In the free 
case by (2.45) £m ax = (ndll - \)kT (which compares with the energy-temperature relation 
E = (nd/2 - l)T(E) in the microcanonical ensemble for the free case) and a -1 = (nd/2 - 
l) -1 ^^,. = (nd/2 - l)(kT) 2 . So as a function of the particle number n in the free case the 
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variance oc~ 1/2 goes like n~ 1/2 when the energy is fixed and like n m when the temperature is 
fixed. We expect this feature to extend to the general case. This would in particular imply 
that at fixed given energy the difference between the microcanonical and the canonical 
description decreases with increasing particle number. 

For periodic and bounded potentials the following asymptotic expansion holds for 
small P 



(2.46) zCfcflS ) = z®;HZ)[i-M + 



p 



2 r 



nV 2 + n(n- l)T 



+ 0(p 3 ) . 



Note that the Laplac e tran sform of the asymptotic expansion ( 2.38 ) gives this asymptotic 
expansion. Relation (2.46) results in the following asymptotic temperature-energy relation 



(2.47) 



nd , 



v 2 -v 2 



+ 0(p 2 ). 



So for fixed small p to leading order the mean energy increases when a potential with 
V > is switched on. Oth erwise it decreases. This behavior is in agreement with the 
corresponding result (2.41) in the microcanonical description. If V = then the mean 
energy increases when such a potential is switched on. 

For large E and T respectively the comparison between the microcanonical and the 
canonical description can be made more explicit. Indeed, relation (2.47) may be inverted 
to give 



1 

T 



(2.48) 



+ e -(P;//<r ) r 2 « 



n + 



nd \ —2 nd^r 



V 



V 2 



which compares with the microcanonical relation ( |2.41[ ). In other words for E and n large 
the temperature defined in the microcanonical ensemble agrees with the one chosen for the 
canonical ensemble, i.e., we have 



(2.49) 



T w r(e(p= \lkT;Hf )\ 



E » e(^=llkT(E;Hf);Hf). 
From (2.47) we obtain the asymptotic expansion for small P 

+ 0(p 3 ). 



knd kn§ 2 



c v (P; = — + 2 



v 2 -v 2 



So if a potential is switched on, the heat capacity increases when p is small. Setting 



(2.50) 



A e -(p;// c ( ;°) 



d 2 > 1/2 



(«)\\l/2 



(kT l c y {T-H7)) 



to be the energy fluctuation we obtain for its asymptotic behavior 

P 2 



(2.51) 



/ nd 1 

Yp 



v 2 -v 



3. The quantum theory 



In analogy to the classical case the microcanonical ensemble for a quantum system 
is usually given as follows. Let {JC, H} be a quantum mechanical system, where IK is a 
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Hilbert space and H a Hamilton operator on "K with a purely discrete spectrum. Moti- 
vated by the classical case on might be tempted to define the the following microcanonical 
quantity 

-j=N (E ; H) = JLkxWE-H)) = tr H (8(£ -//)). 
dE dE 



As is well known (see, e.g., [p5[]) this procedure encounters difficulties for the follow- 
ing reason. Since H has a purely discrete spectrum N(E; H) = tr ^c(&(E - H)) is given 
as the number of eigenvalues (counting multiplicities) of H up to energy E. But then 
dN(E; H)ldE is a sum of 8-functions at the eigenvalues of H with coefficients given by the 
multiplicities of these eigenvalues. In other words dN(E;H)ldE is a generalized function. 
On the other hand if H has continuous spectrum then both trj( (&(E-H)) and U^c(8(E— H)) 
do not make sense. This is in our point of view the main obstacle in making the micro- 
canonical concept useful in the quantum context. The usual way out is to consider "large 
systems", where the spacing of the eigenvalues becomes small and to replace dN(E;H)ldE 
by trjf (©(£ + AE —H)— &(E - H)) with AE being small. This, however, leads to another 
dilemma, since usually there is no natural and intrinsic choice for the size of AE. 

The approach we will propose will give a precise meaning to the notion "large system" 
in the sense that the number of particles n only should be sufficiently large. The one- 
particle operators H = we will consider are supposed to have continuous spectrum. 
More precisely, we will consider one-particle Schrodinger operators acting on the Hilbert 
space L 2 (W I ) and of the form 

H = H + V, 

Here Ho is the free Hamiltonian for a particle of mass m > 

h 2 

H =-—A 
Zm 

where A is the Laplace operator on K , the configuration space for the particle. The poten- 
tials V we have in mind are periodic or random as in ( 2.20| ). 



The resulting n-particle operators with no interaction between the particles will be 
denoted by // (n) . They act on the 71-fold tensor product ® n L 2 (R. d ) of L 2 (M. d ), which is equal 
to L 2 (M" d ), if we consider the particles to be distinguishable or on the subspaces (3„+L 2 (M. d ) 
of symmetric (+) or antisymmetric (-) wave functions, when we allow for statistics. 

Let Xa be the characteristic function of any cube AcM'' and correspondingly let %A" 
denote the characteristic function of A" c M" d . Set Va = %aV. Also let Aa be the Laplace 
operator on A with Dirichlet or Neumann boundary conditions. Let Ha = //o,a + Va with 
//(> A = -h 2 /2m ■ Aa and correspondingly H A ' 7) for the resulting n-particle operator. In the 
examples we have in mind all H^' have a discrete spectrum with no finite accumulation 
points such that N(E;H^) as well as the canonical Boltzmann-Gibbs partition function is 
well defined. Furthermore it can be shown that the limit, called the integrated density of 
states, 



l 

A\- 

exists [ |6l| ] and is independent of the boundary conditions [^6j|. In particular in the free 

nd/2 /^_ m ,E■^n^//2 



(3.1) < n{E\H w ) = lim —N(E-M { " } ) 

- A^~ |A|" A 



(n) 

n-particle case n(£;// ( y") = 0for£<0 



ln) 1 / mE \ (2nmE) n 
(3-2) n(E;H" } ) = -. — T = —, f 

r(^ + i)V2^V h" d r(^- + i) 
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for E > 0. This reflects the well known intuitive observation that each quantum state occu- 
pies the volume h" (h = 2nh) in the phase space R 2 '"'. In other words a comparison with 
the classical case ( Z.22| ) gives the relation h nd n(E; Hq) = nCE;//^"',). So also n(E;H^) is 



smooth for E J= and [nd/2] times differentiable at E = (here [a] denotes integer part of a). 
The smoothness of n(E; Hq), i.e., the order of differentiability, increases with the number 
of particles. 

For random potentials, i.e., when we consider H = = Ho + V a , the integrated density 
of states n(E;H m ) is actually a deterministic quantity, i.e., independent of CO for almost all 
CO. The same is also true for the spectrum spec(H) of H m and in particular for its infimum 
inf spec(//(j)) [^J. This will have the important consequence that all other quantities we 
will introduce using the integrated density of states will also be deterministic. So in order 
to cover both periodic and random potentials simultaneously we will simply write n(E; H) 
for n(E; H m ) and similarly n(E; H (n) ) for the «-particle case. 

Obviously n(E\H (n) ) is a monotone increasing function in E and therefore continuous 
at almost all E. Its derivative, if it exists, is called the density of states for Also by 
construction n(E + c; H (n) + c) = n(E;H {n) ) for any constant c and n(XE; %H {n) ) = n(E; H (n) ) 
for all X > 0. Note also that in general n(E; // (n) ) need not vanish for E < 0. 

If the potential is sign-definite, i.e., ±V > such that ±Va > for all A then by the 
min-max principle (see, e.g., jl^]) we have ^fN(E;H^ ) > ^fN(E;Hq : a) for all A and E 
and n. This gives 

(3.3) Tn(E; H [n) ) > Tn(E; H { " } ) for all E if ± V > 
and more generally 

(3.4) n(E; (H + Vi) (n) ) < n(E; (H + V 2 ) (n) ) for all E if V\ > V 2 . 
Alternatively the integrated density of states may be written in the form 

(3.5) n(£; H w ) = hm ^L- t r Lm „ d) ( XAI ,®(E - ffW)). 
Therefore, for periodic V 

(3.6) n(£; #<">) = tr^^^Ag©^ - //<"»)), 

where Ao is the unit cube in M. d . This relation also continues to hold for random operators, 
when the right hand side is replaced by its average over the random variable. Tha t th e 



integrated density of states is actually deterministic is due to the limit relation in (3.5), 
which is called a self-averaging effect. Below we will show how to obtain n(E; H™) for 
all particle numbers n from the one particle integrated density of states n(E;H). 

One-particle scattering theory may now be used to relate n(E;Ho + V) and n(E ;//(>)• 
More precisely let 8(E;Hq + V\,Hq) denote the phase shift at energy E > for the pair 
(Ho + Va,Ho). In other words, if S(E;Hq + V\,Hq) is the scattering matrix for the pair 
(H + V A , Ho) at energy E then det S(E; H + V A , Ho) = exp(2i8(E;H + V A , H )). For E > 
one has [|| @ 

(3.7) n(E;H)-n(E;H )= lim -j— 8(E;H + V A ,H ) = -§(£;#, #o). 

A^~ 7t|A| 71 

Now -5(E; Hq+Va, Ho)Ik may be replaced by the so-called spectral shift function t,(E; Hq+ 
Va,Ho) of Krein (see, e.g., [Qj) which for negative E equals minus the number of bound 
states below E. We use the spectral shift function to extend the scattering phase to all en- 
ergies in such a way that 5(E; Hq + V\, Ho) = -nt,(E;Ho + V A , Hq). With this extension of 
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8(E; Hq + V A , Ho) to all E relation ( 3.7 ) also extends to all E. In particular, for random po- 
tentials 6(E; H, Ho) is deterministic and n(E; H)=8(E; H, Ho)=0 for £'<min(inf spec(H), 0). 

The trace formula for the spectral shift function leads to the following relation which 
is due to Beth and Uhlenbeck [^], 



(3.8) 



tr,2, 



-P(H +V A ) . 



S / e-V E 8(E;Ho + V A 



H )dE, 



p>0. 



The spectral shift density 8(E; H, Ho) in (3/7) may be interpreted as the phase shift at 
energy E per unit interaction volume. In short we will call it the (total) scattering phase 
density. Also (3.7) extends to the «-particle case giving 



(3.9) n(E;H {n) )-n(E;Hi;') 



1 

: lim — — , 

A^~ jt A " 



8(£; H { o n) + V™,Hk"') = -8(£; H { "\ <') 



A") uW 



(«) zj(»)\ 



Here we have set V A " (x)=% A »(x)V ( -"\x) in analogy to the definition of Va. Also n(E; //'"') = 
S(E;H {n) , H ( Q n) ) = for E < n min(inf spec(tf), 0). 

A precursor to |3~7|) due to Friedel ^ ||] (see also [|, ||, |||) in the context 
of a single impurity is of course well known in solid state physics. To the best of our 
knowledge relation (3.7) and its extension (3.9) to the integrated density of states seems 
to be new. Note that in case of two impurities described by two potentials W\ and W2, 
say of finite range and even non-overlapping, one has 8(E,Ho + W\ + W2, Ho) 5(E; Ho + 
W\ , Ho)+8(E; H0+W2, Hq). Therefore the contribution from several impurities is not simply 
the sum of the contributions from the individual impurities. However, contributions from 
different impurities are asymptotically additive, when the distance between them becomes 
large 



Relations ( p.7[ ) and ( J3 - 8| > are in accordance with the well known chain rule for the 
phase shift 



(3.10) 



8(£; H x , H 2 ) + S(£; H 2 , H 3 ) = 8(£; H U H 3 ), 



valid for any triple of Hamilton ope rators for which the scattering phase for any of the 
pairs exist. In particular, the relation ( 3.1C ) implies 8(£; H\,Ht) = -$(E;H2, Hi). Also the 
equalities 



(3.11) 



8(X£; iHx , XH 2 ) = 8(E + c\H x + c, H 2 + c) = 8(£; H U H 2 ) 



hold for any X > and any real c. Relations (3.11) obviously extend to 8(£; H, Ho). The 



monotonicity property (3.3) is related to the well known relation for the scattering phase 
T8(£; Ho + V A , H ) > for all E if ± V > 0, 



and more generally (compare (3.4)) 

S(£; Ho + Viji. Ho) < 5(£; H + V 2 , A , #o) for all E if Vi > V 2 . 

If V has compact support, then 5(E;Ho + V A , Ho) d oes n ot depend on A for all large A 
and hence in that case n(E; H) = n(E; Hq). Below (see ( 3.37 )) we will see that this implies 
n(E; H M ) = n(E; Hq ) for all n. As mentioned in the introduction it would be interesting to 
obtain information from n-particle scattering theory on n(E;H^ + V*"' + where the 
n-particle interaction is given in terms of two-particle interaction potentials. 

From the discussion so far it is clear that the quantity n(E;H^) is a nice quantum 
candidate to replace the distribution tr^c&(E - H) which replaces the classical distribution 
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N(E; H c {). In particular if n(E;H M ) is differentiable such that necessarily dn(E; H^ldE > 
then we may define 

(3.12) s(E . <H in )) = k ^n(E,p 

dE 

to be the entropy per unit (interaction) volume for a system of n particles. We hasten 
to point out that this definition of the entropy needs a choice of an energy unit since 
dn(E;H^)ldE has the dimension of an inverse energy. In other words, s(E;H^) is only 
well defined up to a constant. One way out is to renormalize the entropy additively by 
choosing a given energy E and to replace s(E;H [n) ) by s(E;H ln) ) - s(E ;H (n) ). Another 
way is to consider a relative entropy (see, e.g., [ f76| | for a general discussion) say for the 
pair (H w ,H ( n) ) as s(E; // (n) , //<,"') = s(E; H (n) ) - s(E;H { q ] ). 

In case n(E;H^) is even twice differentiable w.r.t. E a temperature T(E) = T(E;H M ) 
may be defined by 

1 d d 2 n{E;H {n) )ld 2 E 

(3 - 13) f(E) = dE s(E;E 1 )=k ^kmdE 

which is independent of the choice of the energy unit. Obviously T(E) is positive if the 
second derivative of n(£';// ( " ) ) is positive and below we will show that this is indeed the 

case for the Hamiltonians we consider if we re strict the particle number n to be > 4. 

In analogy to the classical case (see (2A)) the temperature T(E) defined by ( 3.13 ) is 
a monotone increasing function at E if n(E; //'"') is three times differentiable and satisfies 
the inequality 

|,n(E;^)4n(E;^))<(|,„(E : W 2 

In this case the heat capacity c Y (E) defined again by llc v (E) = dT(E)ldE in this quan- 
tum mechanical, microcanonical context is positive. In the opposite case, i.e., with the 
inequality reversed, it will be negative. 

As announced in the previous section the fluctuation of the inverse temperature may 
now be introduced in an analogous way by replacing n(E;H~{) by n(E; H^). For this to 
work we have to assume that the particle number n is so large that n(E; H (n) ) is three times 
differentiable. Below we will show that at least for periodic potentials all derivatives up to 
order 3 are then nonnegative. 

Assume that for given E > inf spec(// ( "') we have strict inequality dn(E; H^)ldE > 0. 
Define the probability measure 



1 d 2 n(E';H w ) , 

dv E (E') = r- -dE' 

v ' dn(E;H^)dE dE' 



on the interval (inf spec(// ( "'), E}. Note the close analogy with the classical case ( 2.33 ). 
Continuing this analogy we introduce the function 

< -J— = —, In -^-n(E';HW) < ^ 
~ 7(E') dE' dE' 2 

With ( • )' E denoting the resulting mean w.r.t. c/Ve we have T(E)~ l = k(1~ l )' E . So again an 
inverse temperature fluctuation is given by 



A'(T- l )(E)=(k 2 (T 2 Y E -T(E) 



-2\ 1/2 



The relation (2.36) involving this fluctuation and the heat capacity again carries over by 
simply replacing n(E;H^) by n(E;H w ). 
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in) 

As a first example we consider the free case (with or without statistics) and for 
which the microcanonical entropy (3.12) becomes negative for small E > and nd > 2 by 
(3.2). As in the classical case by (3.2) and (3.13) the resulting temperature is related to the 
energy by 



(3.14) 



nd 
~2 



1 kT(E). 



For nd large this gives for the energy per particle Eln » dll kT as should be expected. 
Again as in the classical case only when n = d = 1 the temperature so defined becomes 
negative. For nd > 6 the inverse temperature fluctuation is easily calculated to agree with 
the corresponding result (2.37) in the classical case. 

The next relation gives the temperature difference between that of the theory with a 
potential and that of the free theory at the same energy involving the potential only through 
the scattering phase density, 



(3.15) 



T(E;H (n) ) 



T(E;H^) = - 



1 ±S(E; //<">, H<?) -E-&ME; H^, H™) 



k c(nd)E ndl2 - 2 + ^t8(£; H^) 



with c(nd) = Y(ndl2 - l)~ 1 K(ml2nh 2 ) ndl2 . It is a consequence of (3.2) and (3.9). It would be 
interesting to analyze this expression and, in particular, the quantity 8(E; H^- n \ Hq ) when 
E approaches the bottom of the spectrum of // (n) . 

Similarly to our discussion of the classical case we want to argue that at least for pe- 
riodic or random potentials, which in addition are twice differentiable, the heat capacity is 
positive for all large energies. For this we use the large energy asymptotics of the phase 
shift [f7l|,|38|, |65[ ] and on which we shall comment below. In the present context this asymp- 
totics takes the form (assuming in addition that V is twice differentiable in the periodic case 
or that the single site potential Vo is twice differentiable and has support strictly inside the 
unit cube in the random case) we have 

(3.16) _j_8(£; + V$\ //< n) ) = j^n(E; H{?) (ai(n)iT 1 + a 2 (n)E- 2 + 0(E~ 3 )) 



|A| 



with 



fli(n) = 



a 2 (n)-- 



nd 

~Y 

nd 

T 

nd 

T 



nd 
~2 



i ) / v%\xy t <r a x 



nd 1 



2 jnd _ 



V A (x)d d x, 



nd 
~2~ 



— I Y x {xfd d x + 



n(n — 1) 



V A (x)d d x 



What has been proved so far is, e.g., that for periodic or random potentials in one dimension 
there is an upper bound for the phase shift density of the form const ■ E~ m [j47p. In three- 
dimensional case the three-term asymptotics of the integrated density of states is known 



(3.17) 



By the Schwarz inequality 
1 

A\2 



/ Vd.x)d d x 



with equality if and only if V is constant on A. Both terms stay bounded when A 



2 

Let V and V 2 , respectively, denote these limits as A — > °°. They agree with the quantities 
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given by (2.39). In the random case this limit is deterministic, a trivial example of self- 



averaging. 

Combined with ( |3.2| ) and (3.9), and under the assumption that the asymptotic expan- 
sion is preserved in the limit A — > oo, relation (3.16) gives the following asymptotic behav- 
ior of the integrated density of states 



n(E;H^ + V w ) 



nd— 



(3.18) 



E 2 nd I nd 

+ 1 

2 2 2 



nV 2 + n(n- 1)V 2 



+ 0(E~ 3 ) 



Observe that h does not appear explicitly in this asymptotic expansion. Therefore it makes 
sense to compare it with the corresponding classical case. Indeed, when we multiply ( 3.18 ) 
by h" d the resulting derivative w.r.t. E agrees with the derivative of n(E; H^'), This follows 
(3.2) and (2.38) and has the important consequence that to this given order the asymptotic 
energy-temperature relation (2.41) as well as the asymptotic relation (2.42) for the heat 
capacity and the conclusions thereof carry over to this quantum case. Moreover, it is natural 
to conjecture the classical limit relation 



(3.19) 



]imn(E;H in) )h 

no 



nd . 



n(E;H$\ 



For relations in this vein when the classical theory describes ergodic motion, see, e.g., 
[ PH [73|, \7% . In addition, by the example in the previous section o f a cl assical model with 
negative specific heat, with an appropriate control of this limit in (3.19) this could lead to 
a quantum model with negative specific heat for small H. 

We turn to a discussion of the canonical ensemble in the quantum case. In analogy 



to the classical case (see (2.15)) the Boltzmann-Gibbs canonical partition function for H 



(») 



may also be given in terms of N(E;H^ ) as 

(3.20) Z(P; fff) = tr i2(An) ( e -0<) = J e^ E dN(E;H^) = p J e~^ E N(E; Hf)dE 
with mean energy 

d$ Z(p;//f) 

Again by the Schwarz inequality it is a monotone increasing function of the temperature. 

So as in the classical case when A — > °° it therefore makes sense to consider the 
following partition function per unit interaction volume 

(3.21) 

<z(P; H (n) )= hm j^Z^-H^) = J e^ E dn(E; H (n) ) = p J e^ E vv{E;H {n) )dE. 
The mean energy is then given as 

e(|3;//M) = -A i n z( p. ffW) 

(3.22) 



jEe-V E dn(E;H M ) 



1 

-77 + 



z(P;# (n) ) 
P 



Ee-V E n(E; H w )dE 



P z(p;H<")) 

which is not necessarily positive since n(E; may be non-zero for E < 0. Again by the 
Schwarz inequality the mean energy for the canonical distribution is a monotone increasing 
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function of the temperature. For the free case we obtain the familiar relation 



(3.23) 



z(P;#o w ) 



\ nd/2 i 

f m \ 1 

~ h n 



\2%h 2 p ) 



(2nmkT) 



mill 



and hence the energy-temperature relat ion e(p; Hq ) = (nd/2) kT w hich agrees with the one 
in the microcanonical ensemble ( 3.14 ) when n is large. By (3^9) (compare also with the 
relation (3.8)) in the one-particle case, we obtain the relation 



(3.24) z(p;i/ w )-z(p; H { '°) =- [ e^ E d8(E;H M ,H^ } )=^ [ e^ E 8(E;H ( "\H^)dE. 

% J % J 



The relations ( 3.23 ) and ( 3.24 ) give 



e -(P;H w )- e -(P;< J ) = -^ln 



d , z(P; //<">) 



^P z(P;^ n) ) 



4lnfl + f^f /2 I/e-^S(E;^,^) 
dp \ \ m J K J 



2nh 2 p 



1 + 



27^ 2 p sm/: 



e-^dhiE-H^M^) 



(3.25) lj {^-p- l -E\e^ E dh{E-H { "\H^) 



This relation corresponds to the temperature difference < J3 . 1 5| > in the microcanonical de- 
scription and again the right hand side involves the potential only through the scattering 
phase den sity. W e do not know at present under which conditions one may take the limit 
T — » in ( 3.25 ) and what the result should be. As in the microcanonical case this presum- 
ably requires an analysis of the scattering phase density near t he bo ttom of the spectrum of 
//(") which equals n min(inf spec(H), 0) (see the remark after ( 3.32 ) below). 

We note in this context that a theorem of Fumi | [23| | (see also [22|, [jj], |55||) relates 
the shift of the ground state energy due to a single impurity of a system of non-interacting 
fermions to the scattering phase caused by this impurity by zero temperature. We may 
recover Fumi's theorem now in the form of a statement on the shift of the mean energy 
density for non-i nter acting electrons, moving in a periodic or random potential V, from the 
general relation ( 3 . 7 ) in the following way. For given V and given electron density n let the 
Fermi energy Ep = Ef(r\;H) be given as the solution to the equation n/2 = n(£>; H), where 
the factor 1/2 accounts for the spin. Thus E F (f\;H ) = (r\h d T(d/2 + l)/2) 2ld /(2nm) for the 
free theory. Then the ground state energy density Clf is given as 

rE F (r\-H) rE F (n-H) 
(3.26) Q. F (n;H) = 2 / Edn(E;H) = -2 / n(£; H)dE + E F (r\;H)n. 

J inf spec(/f ) J inf spec(//) 

Here again the factor 2 accounts for the spin and the last relation follows by partial inte- 
gration. For the free theory we have Q. F (n;H ) = dl(2(d/2 + l))E F (n; Ho)n. 

We recall that n(E; H) is a monotone increasing function in E but for periodic V it 
is constant outside the spectral bands, so E F may not be uniquely defined if n is such 
that E F (r\;H) lies outside the bands. But by the same observation and the first equality in 
(3.26) Q. F (n,H) is still unique. For random potentials with no gaps in the spectrum the 
deterministic quantity n(£'; H) is strictly increasing, so there is no problem then with the 
definition of E F (n;H). 
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From (3.26) we obtain the following relation for the shift of the mean energy density 
w.r.t. the free theory as 

Q. F (n;H)-a F (r\;H ) = 2 Edn(E;H)-2 / Edn(E;H ) 

iinfspec(H) JQ 

fE F (n;H) fE F (n;H ) 

= -2 n(E;H)dE-2 / n(E;H )dE 

Jinfspec(ff) JO 

+(E F (n-H)-E F (n;Ho))r\ 

2 fE F (n;H ) 
= — / 8(E;H,Ho)dE 

Jt Jmin(inf spec(f/),0) 

in the approximation where E F (n;H) = E F (n;Ho), which should be valid for weak poten- 
tials. Again we note that as in the case of the extension of Friedel's theorem the contribu- 
tion by several impurities is not just the sum of the contributions of the single impurities. 
The correct analysis of the contributions by random impurities was performed in our work 

Returning to our general discussion of the canonical ensemble an entropy may be 
defined in the following way. With p(P;H A n) ) = expi-fiH^IZifcH^) being the density 
matrix (i.e., < p(p;// A ') and tr L 2 (A „)P(p;// A " ) ) = 1) associated to the canonical theory for 
// A at temperature T, let 



< S(p; // A n) ) = -k tr tfCA „) (p(p; flf ) In p(p; ff « ) 

= HnZ(P;<)+— ( Ee-* E dN{E;Hf) 



z(P;#n 

-k I Z(P;#f yV^ln fz(P;fff y'e-^) dN(E;H% } ) 



(3.27) = ~/cp 2 ^^lnZ(P;//r)j = k^^F^H^) 



F(p;// A n) ) = -- In Z(p; //<"') 



be the resulting von Neumann entropy [|57|] . Also 

1 
P 

is the Helmholtz free energy. Recall that S(p;// A n) ) > holds bee ause -x\nx > for 
< x < 1 and all eigenvalues of p(p;// A ) ) lie between and 1. By the second relation in 
( }3.27[ ) the limit 

(3.28) s(p;// (n) )= hm (s(^;H^) - kn In |A|) 

exists and is given as 

s($-H in) ) = ^ In z(p; //'">)+ -Jjj^ J Ee^ E dn(E;H^) 



-k J z(p;// ( " ) r I e- (5£ ln (z(p;i/ ( 'Ve~ P£ ) dn(E;H (n> ) 



(3.29) = -^|(^z(P;^) 

We interpret this as the entropy per unit interaction volume. Since 5(p; // A ) is a monotone 
decreasing function in p so is s(P; //'"') by (3.28). 
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By (3.23) for the free case we have 



nd 



k 1 + In 



fr 2 p 



so in particular s(p;//o"') tends to +°° as (3 



and to — oo as p — > oo. In analogy to 
the microcanonical case ( 3.12| ) the last property is again an undesired feature, which it 
shares with the classical theory (see, e.g., [|76|]). The first property continues to hold in 
the presence of a potential V, which in the present context we will assume to be bounded. 
The way to see this is to consider the difference s (P; H M ) - i(P; H { " ] ). It has an asymptotic 
expansion in p for p small, which is obtained as follows. Using the so called heat kernel 
expansion one can show (see, e.g., [ [25ft and the references given there) that the following 
asymptotic expansion is valid 



1 

|A[»- 



— z(P;nrV 



G^p) 



mill 



V A (x)d d x 



(3.30) 



4 



V 2 A {x)d d x+ n(n 



1) 



V A (x)d d x 



+ 0(p 3 ) 



Here we have assumed that A is a finite union of unit cubes. In addition in the periodic 



case the potential V is assumed to be twice differentiable and in the random case (2.2C) 



the single site potential Vo with support in the unit cube centered at the origin is also twice 
differentiable and vanishes near the b oundary. In case V and Vo have higher order deriva- 
tives the asymptotic expansion ( 3.30 ) also extends to higher orders. Note that apart from 
the pre-factor Planck's constant does not appear. This is a consequence of the assumptions 
just made on V, Vo and on A. In general, powers of h appear combined with derivatives 
of V and V m of the same order within integrals over A. Integration over A therefore gives 
boundary contributions which vanish. 

Assuming again that this asymptotic expansion is preserved in the limit A — > oo, we 
obtain 



(3.31) 



z(P;ff w )=z(P;fl* ,) )(l 



- P 2 — 
■«py + y nV 2 + n(n- 



l)V 2 



+ 0(p 3 ) 



Again up to a fact or h nd this agrees with ( 2.46| ). Using ( 3.21 ) we see t hat (|3 .31 ) is com- 
patible with ( 3.18 ), i.e., the Laplace-Stieltjes transform of ( 3.18 ) gives ( 3.31 ). As a matter 
of fact, the conjecture ma de in [71 1 on the asymptotic behavior of the phase shift (which 
in the present context is ( 3.16 )) was just made to fit in this way with the "heat kernel" 
expansion at large temperatures (small p) of 

1 



-tr, 



For comparison note that in general no asymptotic expansion beyond the famous Weyl term 
for the number of eigenvalues of an elliptic differential operator on a compact manifold like 



N{E;H ( A ) is known. The best result in this direction is due to Hormander J34|]. Although 



related, the conjecture (|3.16[) is of a different nature since instead of the operator it 



involves the operator Hq + V A n } . 

Observe that apart from the overall factor so far h does not appear in the asymptotic 
expansion of z(p; //'"') to this order. This is of course in agreement with our corresponding 
result in the microcanonical setup. Although we shall not need it, we remark that for 
finite range, smooth potentials W in d dimensions the trace of the heat kernel exp{-p(Ho + 
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W)} multiplied by h d and for fixed (3 has a well known asymptotic expansion in (positive) 
powers of h (see, e.g., [^p[|) which has been shown to be a complete asymptotic expansion 
in ||3l[, [72]]. For the asymptotics of the integrated density of states for periodic Schrodinger 
operators we refer to | ]40[ ] . 

Relation (3.31 ) results in the following relation for the entropy 



(3.32) s(P; H { ">) - s(P; //<"') = -jfep ; 



d 1 



z(p;// ( "') 
«/PP"z(P;ff w ) 



•In 



By (3.17) the first term on the right hand side of (3.32) is strictly negative if V 0. We do 
not now about a similar asymptotic expansion in 1/p, which is uniform in A. The problem 
arises because as mentioned above it is difficult to establish the asymptotic behavior of 
(3.25) for small T. Note that for fixed A one has, e.g., 



In Z(p; ff«) 



dim Pi R ,„ . ^ 
-Pe + In dimP + L e -P(£i-e ) 

dimP 



+ 0(e 



). 



where 80 <6] <e2<. . . are the different eigenvalues of Hj"' and P\ the orthogonal projections 
onto the corresponding eigenspaces such that dimP, are the multiplicities of the eigenval- 
ues. For fixed i the difference £,• - £0 tends to zero at least as | A| ud as A — > 00. 
Define 



/(p; H (n) ) = ~ lim (P(p; j n) ) - n|A|) = ~ In z(p; H (n) ) 



Pa 



P 



to be the Helmholtz free energy per interaction volume. From (3.22) and (3.29 ) we obtain 
the relation between the mean energy, the entropy and the free energy, familiar in the usual 
context of the canonical ensemble (compare with the same relation (2.43) in the classical 
case) 

e(P; H {n) ) = rs(p; H w ) + /(p; H {n) ). 



(3.33) 



By (3.30) we have 



/(P;#W)-/(P;#f) = «V-4 



2 



+ 0(p 2 ) 



for small P and where /(P;//q M) ) may be read off (3.23). The mean energy has the asymp- 
totic form 



e(P; //<">) =^P"' 



+ nV- 



np 

Y 



v 2 -v 



+ 0(p 2 ) 



which agrees with the asymptotic behavior (2.47) in the classical canonical description. 
Therefore the inverse asymptotic relation (2.48) carries over and so for large E and n these 
microcanonical and canonical descriptions are approximately equal. In other words in 
analogy to (2.49) we have the approximative relations which are inverse to each other 

T w r(e(P = 1/kT; H {n) )), E w e(p = llkT{E; ff (n) ); H (n) ). 

Also in this canonical setup the heat capacity is defined as 

(3.34) 



c v (p; ffW) = Ag(p ; //(«)) = _L Ae -(p; //W) 2 , 



which is >0 by the d iscussion above and which equals A:nc//2 for all P in the free case as it 
should. As in ( 2.50 ) 

1/2 



(3.35) 



Ae(p;/P'") 



lnz(p;// ( '") 



26 



VADIM KOSTRYKIN AND ROBERT SCHRADER 



is the energy fluctuation. So for small B we obtain 



(3.36) 



c v (P;// (n) ) 



hid knB 2 



V 2 -T 



+ 0(p 3 ) 



and thus for small B the heat capacity increases when an external potential V is switched 
on. By (3.34) relation (3.36) in turn gives the following asymptotic behavior for the energy 
fluctuation 



Ae(B;# w ) 



Ind 1 / B 2 



V 2 - 



-V 2 



+ 0(B 3 ) 



counterpart (2.51) 



which by (3.35) may also be obtained from (3.31) and which agrees with its classical 



In analogy to our analysis in the classical case (see the discussion following (2.44)) 



and by the same arguments (involving the assumption that the heat capacity is sufficiently 
large and positive), also in this quantum context the canonical and the microcanonical 
ensemble both give the same energy-temperature relation and the same heat capacities. In 
particular the varia nce a of the approximating Gauss distribution is just th e squ are of the 
energy fluctuation ( 3.35 ), which in turn is related to the heat capacity by ( 3.34 ). Indeed, 
by observing the necessary differentiability conditions the discussion can be carried over 
verbatim. 

We now want to establish that for given one particle theory with a periodic external 
potential smoothness of n(E; H^') increases with the particle number n and that the deriva- 
tives are all nonnegative. For this the next relation will become relevant. Observe first that 
for all n > 2, 1 < k < n - 1, A and E 



N(E;H™): 



N(E - E'\ H {k) )dN(E'; H ( "~ k) ). 



This follows easily from the identity 

&(E-e-e')= [ &(E-E'-e)8(E'-e')dE'= I &(E -E' -e)d@[E' -e') 



andfrom//i n) = //r(g)I 



ilk) , 



A 



fe(A»-*) + Iz,2(A*) ® h a and Kl\a*)( a ®B) = tr L 2 (A i ) (A)tr i 2 (A „-i ) (B). 
Therefore the integrated density of states for the n-particle theory obeys a similar relation, 
i.e., it may be written as a Riemann-Stieltjes convolution 



(3.37) n(E;H (n) )= / n(E - E'\ H {k) )dn(E'; H {n ~ k) ) ■■ 



■n(-;H {k) ) * s n(-;H ( "~ k) )(E). 

In particular n(E;H^) is an (« - l)-fold iterated Riemann-Stieltjes convolution of the one- 
particle integrated density of states n(E;H^" =l> ) with itself. Here we used the assumption 
that V > -c for some c > 0, such that // (n) > — nc and hence n(E,H (n) ) = for all 
E < -nc. Let E^ esh > -cn be the infimum of all E such that n(E; Z/ (n) ) ^ 0. Since n(£'; H) 
is Holder continuous and nonnegative, the Riemann-Stieltjes convolution of n(E,H) with 
itself exists and is also Holder continuous. This is easily established with techniques used 
in, e.g., Relation ( |3.37| ) implies 

z(B;tf w ) = z(B;H w )z(B; //<"-*>) 

giving in particular z(B; // (n) )=z(B; H) n and the obvious relation /(B; // ( "')=«/(B; H) for the 
free energy. Relations of this type are of course familiar in the ordinary quantum canonical 
ensemble formulation of noninteracting particles. Thus we have, e.g., for the partition 



function given by (3.20) 



Z(B;//f) = Z(B;Z/r*')Z(P;tfD 



(n-k). 
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With reasonable assumptions on n(E,H) concerning the van Hove singularities when V is 
a periodic potential and which hold in d = 1 dimension and for the Kronig- Penney model 
(see, e.g., [ j]], for a survey of spectral properties for periodic potentials in general see [p9|]) 
we will show in the appendix that n(E;H < - n=T> ) is actually continuously differentiable except 
for a discrete set of energies (without finite accumulation points) where, however, the right 
and left derivatives exist. As a consequence n(E; // (,1=3) ) is continuously differentiable at 



all energies. Using (3.37) by complete induction we obtain that n(£; H in) ) is (n - l)/2] 
times continuously differentiable in E with the representation for the /'th derivative (I < 
[(«-l)/2])as 



n (')(E;H w ) 



(3.38) 



n ( '" V - # ( "~ 2 Vn'(£'; # (2) ) 
('-!)(.; * s n '(.;// (2) )(£). 



Also from ( 3.38 ) it follows again by complete induction in / that these derivatives are 
all nonnegative. In particular n(E;H w ) is concave in E for n > 5. This increase of 
smoothness with n is related to the known increase of smoothness with the dimension 
d of space. Indeed write Hq to denote the dependence of the free Hamiltonian on the 
space dimension. Then for given external one particle potential V we have (Hq + V) (n) = 
H$ d + V ( "\ where as in the classical case the potential V (n) is given as V ( "'(x) = V(x,) for 
x = (x\,X2, ■ ■ ■ , x n ) e W d and with jc,- e W 1 . Note that periodicity extends in the sense that 
e.g. V {n) (x + j) = V {n) (x) holds for any j e Z nd and x e W d , whenever V(y + i) = V(y) for all 
i e Z f/ and y e R rf . 

To sum up by our preceding discussion for one particle Schrodinger operators with 
periodic potentials we therefore may define a temperature for the microcanonical ensemble 
for all H (n) , n > 5. Since the third derivative of n(£';// ( "'),n > 7 is non-negative it 
can not vanish identically on any interval [E^ esh ,E]. To see this, assume the contrary. 
Then by integrating 3 times we see that n(£';// ( " ) ) also would vanish identically on this 
interval, which is not possible. So again by integrating, the first and second derivatives 
of n{E\H (n) ),n > 5 are strictly positive for any E > E^ esh . So we conclude that the 
microcanonical temperature for H (n) ,« > 5 is actually positive, finite and a continuous 
function of E. Similarly, since the fourth derivative of the density of states for H (n) , n > 9 
is non-negative, the second derivative is concave in E. However, this is not sufficient to 
conclude that 1IT is a concave function of E whenever n > 9 and presently we do not 
know sufficient conditions on the potential V which ensure this. 

We expect similar results to be valid for the case of stochastic potentials, where it is 
known that the integrated density of states is Lipschitz continuous, i.e., Holder continuous 
of index 1 (see [||]), so one has smoother properties than in the periodic case. For easier, 
recent proofs that they are Holder continuous of any index smaller than 1, see [ fl3| , [jgl ]. 
Unfortunately, in general the Riemann-Stieltjes convolution does not improve Holder con- 
tinuity or any other similar kind of regularity (see, e.g., [24, 17|]). It would be interesting to 
find additional properties of the integrated density of states which allows one to conclude 
that their Riemann-Stieltjes convolution improves regularity. 

So far we have not taken Bose-Einstein or Fermi-Dirac statistics into account. To 
do this observe first that the permutation group S„ of n elements has a canonical unitary 
representation 7t — > U(k), k e S„ on (g)„Z, 2 (R d ). A similar observation is valid when R d 
is replaced by A. The orthogonal projections onto the subspaces of symmetric (+) and 
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antisymmetric (-) wave-functions ® n ,±L (K ) are given as 

(3.39) P { " ] = 4- E (sign n) ( ^ m U {%) 

respectively, i.e., ® n ,+L 2 (R d ) = P| n) <g>„ L 2 (R d ). This gives ti^ i±LHm (A) = tr^^P^A) 
and similarly with M. d being replaced by A. Then one may show that the resulting density 
of states are given by 

(3.40) n ± {E;H (n) ) = Urn v^tx L z m (if >0(£ - #| n) )) = ~n(E; // ( " , ) > 



i.e., when inserting ( 3.39 ) into the middle term in ( 3.4C ) the contributions from K^id vanish 



in the limit. In the periodic case and in analogy to (3.6) one also has 



n±(£; H {n) ) = tr LW %K P ± ®( E " //( '") 



In particular the resulting temperature in the microcanonical ensemble does not depend on 
the statistics. Similarly the partition function in the canonical ensemble is modified by a 
factor 1/n!. 

To conclude this section we briefly describe the associated grand canonical ensemble 
for both statistics using the integrated density of states. In a finite volume A the partition 
function for the grand canonical ensemble is given as 



«=i v 

with /j being the chemical potential. Standard arguments therefore give 

lnZ ± (p>;ff) = ^(lT^I^ / (iTe-^Y'dNiE-M^l 

Hence z+(p, J u;//) may be defined in terms of the one-particle integrated density of states 
by 

lnz ± ($,n;H)= lim -j- lnZ+ft.^H) = [ ( lTc* ,; ) Tl ^ln(£;^/ ( ° =l, ). 
a^~ |A| J V / 

The usual relations for the mean energy and the mean particle number may now be obtained 
easily. 

Appendix A. 

In the the first part of the appendix we will establish increasing regularity for the quan- 
tum integrated density of states with increasing particle number in the case of a periodic 
potential. In the second and last part we will show in an example how randomness may 
serve to increase the smoothness of the classical integrated density of states already in the 
one particle case. 

For periodic potentials V, where bands appear, it is well known that the one-particle 
integrated density of states n(£'; // (n=1) ) is constant for E inside the gaps, smooth inside 
the bands and behaves like \E - Ej + \ 112 at an upper or lower edge E;_ of a band, so in 
particular it is Holder continuous of index 1/2 everywhere. More precisely we make the 
following assumption, which has been proven rigorously in d=l dimension, see [63]. Write 
the bands (finitely or infinitely many) as the closed, pairwise disjoint intervals [£', _,£,■ +]. 
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There are smooth functions n,- +(E; H), E > small, with «; +(E = 0; H) = such that for E 
near E,+ 



(A.l) n(E;H)- 



\E - E it± \ m n h± (T(E - £,+);//) + c,-,± if T(E - E ii± ) > 
Ci,± if ±(£-£,+)>0 



where c, + = n(£',,±; H). In particular n(£; H) is Holder continuous of index 1/2. Also away 
from the ends Ej + of the bands the function n(E; H) is supposed to be smooth. We note that 
for d > 1 the Bethe-Sommerfeld conjecture is actually true: As established by Skriganov 
and Karpeshina, there are only finitely many gaps (see, e.g., [^9|] for an extensive list of 
references and [fffi]). 

With these assumptions we will show that vv(E;H^ n=T> ) is continuously differentiable 
in E from the right and from the left. Also the right and left derivatives agree except for 
the discrete set of energies of the form E = E t + + Ei+. In other words the derivative of 
n(£'; // ( " =2) ) is possibly discontinuous at these energies. Furthermore we will show that 
n(£'; // <B=3) ) is differentiable in E for all E with a derivative which is Holder continuous of 
index 1/2. 

For the proof let < %, +(£) < 1 be smooth functions which are equal to one near 



Ei+ and where (A.l) holds. They may be chosen to have non-overlapping support, i.e., 
any product of two different X/,±(£)' s is zero. Set < %(E) = L,,±Xi,±(£) < 1 an d write 
ni(£; H)=(l-%(E))n(E\ H), n 2 (E;H)=%(E)n(E\ H) such that n(E; H)=n 1 (E;H)+n 2 (E;H). 
Then ni (E;H) is smooth everywhere and xi2(E;H) vanishes outside small intervals around 
the points E, ■ +. With this decomposition of n(E; H) we obtain 

(A.2) 

n(£; H 2 ) = m(-; H) * s m(-; H)(E) + 2m(-; H) * s n 2 (-,H)(E) + n 2 (-; H) * s n 2 (-; H)(E). 



The first two terms in ( |A.2j ) are easily seen to be continuously differentiable w.r.t. E, 
so we only have to show that the third term is continuously differentiable. It suffices to 
consider any term of the form 



(A.3) J XidE ~ E')n(E - E'; H)d{%j,±(E')n(E'; H)), 

since for given E there are only finitely many Ej±s with Ej+ < E. Now we use the 



assumption (A.l) and obtain 

d(x,, ± (EME;H)) = &(t(E-E l± )) {\E - E l , ± \ m f lxl {E)+\E - E L± \- m f l± , 2 {E)) dE 
+fi,±,3(E)dE 

I = i or j, where the //,±,*'s, k = 1, 2, 3 are smoo th w ith support contained in the support of 



X/,±. We insert this into the formal derivative of ( |A.3| ) and obtain 9 terms all of which should 
be finite integrals. The potentially most dangerous one and on which we shall concentrate 
is 

(A.4) 

J &(T(E -E'- E i± ))®(T(E' - E l± ))\E -E'- E^^E 1 - E h± \- m f i± , 2 {E - E')f jX2 {E')dE' . 

For given i, ± and j, ± this integral is finite and continuous in E when E stays away from 
Ej + + Ej+, so let E = Ej + + Ej + + e with e small. Now we make the variable transformation 



E' = Ej + + e'. Then ( A4) may be written as 



(A.5) j 0(t(£ - e'))0(Te') |e - e' |" 1/2 |e' |" 112 gi (e - e')g2(e')de' , 
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where gi(e— 6')=/;,±,2(£ — £'+#;,±) and g2(£ ! )=fj,±,2(£'+Ej±) are smooth and vanish outside 
a neighborhood of the origin. Also the choice of the sign in =p in the first Heaviside step 
function refers to i, ± while the choice of the sign in =F in the second refers to j, ±. For e^O 
(A.5) is continuous in e. Also both limits £ J, and £ f exist but are in general different. 
In fact, consider the case +,+ of the signs, i.e., the choice i, - and j, —. Then (|A.5[ ) vanishes 
for e < and for e > we make the variable transformation e' = ex such that ( A.5 ) takes the 
form 

\l-x)- m x- m gl (e(l -x))g 2 {ex)dx, 
which is continuous in e > and has a finite limit equal to 

*i(0)* 2 (0) f\l-x)- m x- m dx 
Jo 

when e j. 0. The other three cases are discussed similarly. This concludes the proof of our 
claim on the differentiability of n(E; H^" =2> ). To see that n(E; i/ ( "" 3) ) has a derivative which 
is Holder continuous of index 1/2, we use the representation ( 3.38 ) with n = 3 and I = 1 
for the derivative of n(E; Z/ ( " =3) ). Since we just established that n'(E; // ( " =2) ) is continuous 
in E except for a discrete set of discontinuities (without finite accumulation points) some 
easy arguments establish this last claim. 
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